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FORTRAN PROGRAM FOR COMPUTING THE PRINCIPAL MOMENTS 


OF INERTIA OF A RIGID MOLECULE 
By Janet G. Ehlers and Glenn R. Cowgill 
Levis Research Center 


SUMMARY* 

A program to compute the principal moments of inertia of a rigid molecule 
was written in the FORTRAN language. This program eliminates considerations of 
symmetry and avoids hand calculations by reducing the input to structural param- 
eters , that is, bond lengths and angles. The assignment of the structural param- 
eters follows geometric models. In its present form the program is written for 
molecules with not more than 24 atoms. 


INTRODUCTION 

Many structurally complicated molecular species are becoming of interest at 
shock- tube and rocket- combustion temperatures. These include polyatomic species, 
such as HBOg and BFC1, which have no axes of symmetry, and polymers, such as 
(NaOH )2 and (HB02)3* Obtaining thermodynamic functions for these species re- 
quires that the principal moments of inertia or the product of the principal mo- 
ments be known. The calculation, however, becomes cumbersome for large molecules 
or for molecules in which the center of mass and the directions of the principal 
axes through the center of mass are not immediately apparent. In order to elim- 
inate considerations of symmetry and to avoid hand calculations, the method of 
reference 1 together with the authors* technique to minimize the input was pro- 
gramed in the FORTRAN language. 

The program described herein was used to compute the moments of inertia re- 
ported in reference 2. In its present form the program is written for molecules 
with not more than 24 atoms; however, extension to consider molecules with more 
than 24 atoms could be easily accomplished. Any set of atomic weights with the 
corresponding value of Avogadro’s number may be selected as input to the program. 

The program decks are available from the authors upon request in either 
FORTRAN II or FORTRAN IV. Also Included for the user’s convenience are data 
cards for the atomic weights based on the chemical scale of natural oxygen 
(0 = 16.0000) and the 1961 Table of Atomic Weights (see ref. 3) based on the 
exact weight of 12 for carbon 12. 



SYMBOLS 


A 

^-xx’^y-y >-*-zz 

^xy’^yz’^xz’ 

T T T 
J -yx-’ x zy- ,x zx 

M 

% 

r 

r i,k 

X,Y 

a 

A 

CT i 

<Pi,9 k 


symmetric matrix 
moments of inertia 
products of inertia 


i 

mass of the i^* 1 atom 
bond length 

bond length between the i^* 1 and k^* 1 atoms 
atoms in pyramidal or tetrahedral molecules 
Cartesian coordinates of the i^k atom 
angle 

branch angle with vertex at the i^* 1 atom 
eigenvalue 

out- of -plane-atom angle with vertex at the i^* 1 atom 
chain angles with vertex at the i^* 1 and k^* 1 atoms 


METHOD OF CALCULATION 

The principal moments of inertia are the eigenvalues of the real symmetric 
matrix A 



The diagonal elements are the moments of inertia about the center of mass; the 
off-diagonal elements are the products of inertia about the center of mass. 

The matrix elements are calculated according to the equations given in ref- 
erence 1. These equations, where mq is the mass of the i^* 1 atom, xq, yq, zq 


2 



are the Cartesian coordinates of this atom, and M =* ^ m i> are reproduced for 
convenience: 1 


L xx 


L yy 


-zz 


= E®i(yi + 4 ) - M (E“±y i) 2 " s 

= + z i) - i (Svi) 2 - i (E^n)' 

E“i( x i + yf) - i (E% x i) 2 - M (E“±yi)' 

= - 2 m i x iyi + i (E^i) (E m iyi) 

- E Wl (E% x i)(E%Zi) 

- E^yi 2 ! + i (E^/i) (E“i z i) 


T = T 
X xy X yx 


-^xz “ *zx 


%z “ J zy ~ 


Any convenient set of orthogonal axes may be chosen to describe a molecule. 
The equations for the matrix elements give the moments of inertia and the prod- 
ucts of inertia about a set of axes parallel to the assigned axes and with origin 
at the center of mass* 

A convenient value in the calculation of thermodynamic properties is the 
product of the principal moments of inertia. This value is simply the determi- 
nant of the matrix A. 

The three principal moments of inertia can be found individually as the 
roots A of the secular equation 


■xx ~ ^ 

I xy 

-J-XZ 

■yx 

*yy " A 

*yz 

•zx 

I zy 

3*zz 


= 0 


The derivation of this equation is discussed in texts on classical mechan- 
ics, such as reference 4. 
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An attempt to solve directly for the roots of the preceding cubic equation 
resulted in the appearance of false imaginary solutions , which were caused by in- 
herent machine inaccuracies and rounding errors. Of several other possible meth- 
ods, the iterative Jacobi method described in reference 5 was chosen to determine 
the principal moments inasmuch as this method does not depend on the existence of 
three unequal real roots. In the solution is found an orthogonal matrix, that 
is, an eigenvector matrix, which transforms the original symmetric matrix A 
into a diagonal matrix. This is accomplished through the annihilation of se- 
lected off-diagonal elements by elementary orthogonal transformations. The ele- 
ments of the resulting diagonal matrix are the principal moments of inertia. 

Because of rounding errors in the computation, the principal moments of in- 
ertia are accurate to approximately six significant figures. This degree of ac- 
curacy was determined by repeating computations for a molecule after a change in 
its orientation relative to the assigned axes. This accuracy, however, is ade- 
quate in the light of the reliability of experimental structural parameter data. 


DISCUSSION OP THE PROGRAM FOR VARIOUS TYPES OF MOLECULES 

Figure 1 is a flow diagram of the main program. The block numbers corre- 
spond to the statement numbers of the code listing. A flow diagram of the sub- 
routine Jacobi is given in reference 5. The computation of the principal moments 
of inertia is simplified by equations for the Cartesian coordinates of atoms in 
terms of the appropriate bond lengths and angles. A molecule is classified ac- 
cording to its geometric type; planar, pyramidal, tetrahedral, other nonplanar, 
or general. The input for each geometric type follows a specified model. Al- 
though this input is reduced to bond lengths and angles, their order of assign- 
ment is thus not completely arbitrary. 

Each of the six geometric types is illustrated in figure 2 with the equa- 
tions for the atom coordinates given in table I. Sample computations for each- 
type of input are given in table II. These computations can be used to check the 
program deck. 

Table III is a listing of the code preceded by an explanation of the input 
cards, and followed by a listing of the input cards used for the computations in 
table II. The atomic weights and the corresponding Avogadro’s number are based 
on the chemical scale of natural oxygen. 


Type 1 - Planar 

The input for planar molecules can be reduced to bond lengths and angles 
when the molecule is considered to be a chain of atoms with other atoms branching 
from the chain. Planar molecules without and with branching are illustrated in 
figures 2(a-l) and (a-2), respectively. Figure 2(a-l) represents a molecule 
such as HBO 2 , which has a chain of four atoms; figure 2 (a-2) represents a mole- 
cule such as (HB 02 ) 3 > which can be visualized as a chain of six atoms with three 
two-atom branches attached to the chain. The branch atoms are specified accord- 
ing to the chain atom from which they branch, that is, atom 10 is the first 
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I 


branch atom attached to chain atom 1, and atom 30 is the first branch atom at- 
tached to chain atom 3. 

The Cartesian coordinates of the atoms are computed by the program from the 
equations in table 1(a). Atom 1 is considered to be at the origin of a set of 
orthogonal axes and bond r 2 3 to be parallel to the x-axis. The general equa- 
tions given in table I (a) will permit the extension of the program to accommodate 
more chain and branch atoms. Presently, the program can treat a nine-atom chain 
with two-atom branches from any of the first eight chain atoms. 

The necessary input data are the appropriate bond lengths and angles accord- 
ing to the following assignment: 

(1) The adjacent bond lengths are specified in order along the chain start- 
ing from atom 1 with bond r^_ 2 . 

(2) The chain angles are specified in order along the chain beginning with 
angle cp 2 between the first and second bonds. The assigned angles are always on 
the same side, that is, the "inside, 11 of the chain. 

(3) For any particular branch, the bond lengths and branch angles are speci- 
fied in order from the chain; all angles are taken on the same side of the 
branch, as indicated in figure 2(a-2). 

Although (HB0 2 )3 is a closed cyclic structure, it is illustrated in figure 
2(a-2) as an open chain with the bond length rg and angle cpg omitted. 

These parameters are unnecessary because the preceding bond lengths and angles 
have specified all the atom positions. 

Other chain assignments are also possible for for example, an 

eight- atom chain with origin at the illustrated atom 11; the two branches would 
then be attached to the fifth and seventh chain atoms. 

The equations for the atom coordinates do not permit branching from the last 
chain atom. If this error is made in the branch assignment, the computation 
stops and a suggested assignment is printed. 

The input for the simplest chain, a diatomic molecule, consists only of the 
bond length. 

Sample computations for planar molecules without and with branching are 
given in table II (a). 


Type 2 - Pyramidal XI 3 

Several molecules such as NH3, FH3, PCI3, and FF3 are of the pyramidal XI 3 

geometric type. The input for this type can be simplified to the X-I bond 
length and the Y-X-I angle. 
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A three-dimensional pyramidal molecule, oriented relative to a convenient 
set of orthogonal axes, is illustrated in figure 2(b). The figure also includes 
diagrams of the atom projections on the xy-, xz-, and yz-planes. The X atom 
is at the origin of the orthogonal axes and the Yj_ and the Yg atoms lie in 

the xz-plane. The z-axis bisects the Y^-X-Yg angle. 

The equations for the Cartesian coordinates of the atoms illustrated in fig- 
ure 2(b) are given in table l(b). Because of the selection of the coordinate 
system, these equations are simple functions of the X-Y bond length r and the 
Y-X-Y angle a,. 

A sample computation for a pyramidal XY 3 molecule is given in table Il(b). 


Type 3 - Tetrahedral XY4. 

Molecules such as CCI4, SiCl^, SiF^, and SiH^ are of the tetrahedral XY 4 
geometric type. The only input necessary for this type is the X-Y bond length. 

A three-dimensional tetrahedral molecule, relative to a convenient set of 
orthogonal axes, is illustrated in figure 2 (c). The figure also includes dia- 
grams of the atom projections on the xy-, xz-, and yz-planes. The *1* Y 2> and 

Y 3 atoms lie in the xz-plane with Y 2 on the z-axis, and the X and Y4 atoms 
lie on the y-axis. 

The equations for the Cartesian coordinates of the atoms illustrated in fig- 
ure 2(c) are given in table 1 (c). Because of the selection of the coordinate 
system, these equations are simple functions of the X-Y bond length r. 

Sample computations for the tetrahedral XY4. molecules are given in table 
11 (c). The input for the P4 computation involves a fictitious X atom, which is 
discussed in the section, SPECIAL TREATMENTS OP MOLECULES. 


Type 4 - Other Nonplanar 

Certain molecules have one or two atoms that are not in the same plane as 
the other atoms. Examples of these molecules are (NaOH) 2 , S 2 C1 2 , and SF^. Their 

out-of- plane atoms are the two H atoms, one Cl atom, and two F atoms, respec- 
tively. Three-dimensional diagrams of these molecules, oriented relative to a 
convenient set of orthogonal axes, are sketched in figure 3. The input is speci- 
fied by considering the out-of -plane atoms to be special branches from a planar 
chain. In addition to the bond lengths and angles required for the planar chain 
input, the out-of-plane-atom (hereinafter called O-P-A) bond length and angle are 
necessary. 

In figure 2(d) a three-dimensional nonplanar molecule is illustrated with 
the atoms oriented relative to a convenient set of orthogonal axes. The figure 
also includes diagrams of the atom projections on the xy- and xz-planes. The 
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out-of-plane atom is specified according to the chain atom from which it 
branches; that is, out-of -plane atom 200 is attached to chain atom 2. A single- 
atom 0-P-A branch is permitted only from chain atom 2 or 3. This restriction, 
however, is not severe, and many nonplanar molecules, including all those of in- 
terest in reference 2 , can be described in this manner. 

The Cartesian coordinates of the out-of-plane atoms illustrated in figure 
2(d) are computed from the equations in table l(d). The 0-P-A bond length r^ k 

is the distance between the out-of-plane atom and its adjacent chain atom. The 
0-P-A angle cr^ is formed by the out-of-plane atom, its adjacent chain atom, and 

one other chain atom, provided that these three atoms form a plane perpendicular 
to the plane of the chain. The intersection of these two planes is parallel to 
the x-axis and contains the line segment corresponding to bond length r 2 3 - A 

plus or a minus sign indicates that the out-of -plane atom is above or below the 
plane of the chain. 

Sample computations for other nonplanar molecules are given in table Il(d). 
Special discussions of the input for the (NaOH )2 and the SF 4 computations are re- 
served for a later section, SPECIAL TREATMENTS OP MOLECULES. 

[types 5 and 6 - General; Coordinates Read In 

In these two types the atom coordinates may be read in either in the form 
U + V cos a (type 5) or directly as x,y,z (type 6 ). When the orthogonal axes 
are judiciously chosen, the Cartesian coordinates of the atoms are usually simple 
functions of the bond lengths and angles, types 5 and 6 are useful in the com- 
putation of principal moments of inertia for molecules not included in the other 
types and as checks on the other computations. 

Examples of type 5 and type 6 computations are given in tables 11(e) and 
(f), respectively. The molecules, oriented relative to the assigned coordinate 
system, are illustrated in figures 2 (e) and (f). 


SPECIAL TREATMENTS OP MOLECULES 

A facility in visualization of structure is helpful in choosing the simplest 
way of treating a molecule. This is particularly the case with Type 4 - other 
nonplanar molecules. Occasionally the intuitive order of atoms along a chain is 
not the most convenient. For instance, the order of the chain atoms for (NaOH) 2 ^ 
illustrated in figure 3(a), was taken to be Ra 0 0 in order to accommodate the 
H atoms as 0-P-A branches. 

Certain molecules can be made to fit into a desired geometric- type classi- 
fication by the invention of fictitious atoms of zero mass. The molecule SF 4 , 

shown in figure 3(c), was assigned a fictitious third chain atom in order to per- 
mit a planar branch and two 0-P-A branches from the second chain atom. Thus the 
0-P-A plane is perpendicular to the plane of the chain. 
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Another example of a molecule for which a fictitious atom was invented is 
P4. Ordinarily a Type 2 - pyramidal XI 3, this molecule can be treated as a 

Type 3 - tetrahedral XY4. In this way, the available experimental data are more 

easily accommodated. The X atom is fictitious and the four Y atoms corre- 
spond to the P atoms. The required bond length input is the fictitious XY bond 
length, which is readily calculated from the experimental data. Since a ficti- 
tious atom has zero mass, it does not affect the principal moments of inertia. 


Lewis Research Center 

Rational Aeronautics and Space Administration 
Cleveland, Ohio, Rovember 13, 1963 
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TABLE I. - EQUATIONS FOR THE CARTESIAN COORDINATES OF THE ATOMS 
(a) Type 1 - planar 


Atom 

x-coordinate, 

angstroms 

y~ coordinate, 
angstroms 

1 

0 

0 

2 

-1*1 2 cos q>2 

r l,2 sln f>2 

3 

X 2 + r 2,3 

y 2 

4 

x 3 - r 3<4 cos <p 3 

y 3 - r 3j4 sin <p 3 

5 

x 4 + r 4,5 °°3(?3 + 9 4 ) 

y 4 + r 4,5 sln (<P 3 + 9 4 ) 


/l-l \ 

/l-l \ 

a i 

x i-l + (“ 1 ) 1 “ lr 1 -l,i cos ( Yj 

y l-l + (- 1 ) 1 ~ lr i-l,l 3lr H 2 9k) 


\k=3 / 

\k=3 / 

10 

-'1,10 003 ( e i - *P2 ) 

_r l,10 aln ( 0 l - ^2) 

11 

x io + r io,ii 003 ( e i + e 10 - <p 2 ) 

y 10 + r 10,ll 3ln ( 0 l + 6 10 - 9 2 } 

20 

x 2 + r 2,20 003 9 2 

y 2 + r 2,20 3ln 9 2 

21 

x 20 ~ r 20,21 ooa (®2 + ®20) 

y 20 - r 20,21 sln (®2 + 0 2O) 

30 

x 3 " r 3 ,30 cos(9 3 + 9 3 ) 

y 3 “ P 3,30 Sln ( e 3 + 

31 

x 30 + r 30,31 coa ( e 3 + e 30 + 9 3 ) 

y 30 + r 30,31 sln ( 0 3 + S 30 + 93) 

40 

x 4 + r 4,40 cos (®4 + <P3 + 9 4 ) 

y 4 + r 4,40 Sin ( e 4 + ’3 + *4) 

41 

x 40 “ r 40,41 003 ( e 4 + s 40 + 93 + 94) 

y 40 - r 40,41 aln < 0 4 + ®40 + * S + ^ 

50 

x 5 " r 5,50 003(65 + fs + 9 4 + 9 5 ) 

y 5 “ r 5,50 sln ( e 5 + ¥3 + 9 4 + 95) 

51 

x 50 + r 50,51 cos ( 0 5 + 650 + 9 s + 94 + <P 5 ) 

y 50 + r 50,51 aln ( e 5 + e 50 + 1>3 + 94 + 95) 

b 10i 

x i + ( - 1)ir i,ioi 003 (®i + E 3 <p k j 

y i + (-^S.lOi aln (®l + £ 9 k j 

b 101+1 

x 10i + ( _1 ) r 101, 101+1 003 (®i + ®10i + k ^ 3 9 k ^ 

y 101 + ( _1 ) i lr 10i, 101+1 3in (®1 + ®ioi + 9 k j 


a General equations for extension of chain. 
^General equations for additional branches. 



TABLE I. - Concluded. EQUATIONS FOR THE CARTESIAN 
COORDINATES OF THE ATOMS 
(b) Type 2 - pyramidal XYj 



(c) Type 3 - tetrahedral XY 4 


X 

0 

r/3 

0 

Y 1 

r 

r 3 

0 


y 2 

0 

0 


*3 

-r n/6 
3 

0 


Y 4 

0 

r _l 

0 


(d) Type 4 - other nonplanar 


200 

x 2 + r 2,200 003 0 2 

y 2 

( ± ) r 2,200 8ln ct 2 

300 

X 3 " r 3,300 008 ff 3 

y 3 

^^SOO 8ln ct 3 


a k = 1 - 3 sin 2 % + cos 2 . 
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TABLE II 


SAMPLE COMPUTATION OUTPUTS FOR PRINCIPAL MOMENTS OF INERTIA 


(a) Type 1 
(a-1 ) KB0 2 


PLANAR 

WITH 4 

CHAIN 

ATOMS AND 0 

BRANCHES. 

WHITE JCP V32 

P408 FEB 1960. 



NO. 

CHAIN 

ATOM 

BOND LENGTH 

ANGLE 

MOLECULAR WEIGHT 

X 

Y 


1 

H 


1.0000 


1.00799999 

0 

0 

0 

2 

0 


1.3400 

120.00 

16.0000000 

0.49999996 

0.86602541 

0 

3 

B 


1.2000 

180.00 

10.8199999 

1.83999993 

0.86602541 

0 

4 

0 




16.0000000 

3.03999990 

0.86602537 

0 


THE GENERATED SYMMETRIC MATRIX IS ...... 

I 

I 0.1226276 -0.2531333 0 

I 

I 

I -0.2531333 9.0980780 0 


THE EIGENVALUES, LAMBDA, OF THE MATRIX ARE THE 
PRINCIPAL MOMENTS OF INERTIA (G-SQ CM) 


LAMBDA! I 1*1. 0E*39 * 0.11549414 
L AMBDA ( 2 ) * 1 . 0E + 39 * 9.10521126 
LAMBDA! 31* 1.0E+39 * 9.22070527 


0 


0 


9.2207053 


PRODUCT LAMBDAS*!. OE+l 17 * 9.69648027 


THE GENERATED MATRICES FROM THE JACOBI METHOD ARE AS FOLLOWS, 


THE DIAGONAL MATRIX (Dl 

0.11549414 0 

0 9.10521126 

0 0 


0 

0 

9.22070527 


THE 

EIGENVECTOR 

MATRIX 

(S) 


I 

I 

0.99960316 


-0. 28 1 692 7 1 E-0 1 

0 

I 

0.281 6927 1 E- 

■01 

0.99960316 

0 

I 

I 

0 


-0 

1.00000000 


THE MATRIX CS)*(D)*(S TRANSPOSE), WHICH 
SHOULO EQUAL THE ORIGINAL MATRIX IS, 


0. 12262755 
-0.25313329 
0 


-0.25313329 

9.09807754 

0 


0 

0 

9.22070527 
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TABLE II. - Continued. SAMPLE COMPUTATION OUTPUTS FOR PRINCIPAL MOMENTS OF INERTIA 

(a) Concluded. Type 1 
(a-2) (HB0 2 ) 3 


PLANAR WITH 6 CHAIN ATOMS AND 3 BRANCHES. 


WHITE JCP V32 P488 FEB. 1960 


6-0 1.36A, 0-H 1 . A, 


NO. CHAIN ATOM 

1 B 

2 0 

3 B 

4 0 

5 B 

6 0 


BOND LENGTH 

1.3600 
1.3600 
1.3600 
1.3600 
1. 3600 


120.00 

120.00 

120.00 

120.00 


MOLECULAR WEIGHT 
10.8199999 
16.0000000 
10.8199999 
16.0000000 
10.8199999 
16.0000000 


0.67999993 

2.03999990 

2.71999982 

2.03999972 

0.67999974 


1.17779455 
1. 17779455 
-0 

-1.17779449 
-1. 17779440 


BRANCH NO. 1,- ATTACHEO TO CHAIN ATOM 8 


NO. 

7 

8 


ATOM 

0 


BOND LENGTH 
1.3600 
1.0000 


ANGLE 

120.00 

240.00 


NUMBER 1 


16.0000000 

1.00799999 


-1.35999998 

-1.86000004 


-0.86602537 


BRANCH NO- 2, ATTACHED TO CHAIN ATOM B 


NO. 

9 

10 


ATOM 

0 


BOND LENGTH ANGLE 

1.3600 120.00 

1.0000 240.00 


16.0000000 

1.00799999 


2.71999997 

2.22000015 


2.35558903 

3.22161448 


BRANCH 

NO. 3, 

ATTACHEO 

TO CHAIN 

ATOM B 

, NUMBER 5 




NO. 

ATOM 

BONO 

LENGTH 

ANGLE 





11 

0 

1 . 

.3600 

120.00 

16.0000000 

2.71999946 

-2.35558909 

0 

12 

H 

1 . 

.0000 

240.00 

1.00799999 

3.71999943 

-2.35558921 

0 


THE GENERATED SYMMETRIC MATRIX IS 


THE EIGENVALUES, LAMBDA, DF THE MATRIX ARE THE 
PRINCIPAL MOMENTS DF INERTIA (G-SQ CM) 


44.624305 -0 .370234 3E-05 0 

-0 . 3702343E-05 44.624300 0 


LAMBDA! 1 ) * 1 . 0E + 39 
LAMBDA 121*1. OE + 39 
LAMBDA! 31* 1 .OE+39 


44.6243057 

44.6242971 

89.2486057 


0 


0 


89.248606 


PRODUCT LAMBDAS*!. 0E + 1 17 * 177723.270 


THE GENERATED MATRICES FROM THE JACOBI METHOD ARE AS FOLLOWS, 



THE DIAGONAL 

MATRIX (D) 



THE 

EIGENVECTOR MATRIX 

IS) 



44.6243057 

0 

0 

I 

I 

I 

I 

0.86644093 

0.49927959 

0 


0 

44.6242971 

0 

I 


-0.49927959 

0.86644093 

0 


0 

0 

89.2486057 

I 

I 

0 

0 

1.00000000 


THE MATRIX !S)*!D!*(S TRANSPOSE), WHICH 
SHOULD EQUAL THE ORIGINAL MATRIX IS, 


44.6243029 

-0.35762787E-05 

0 

■0.40531 158E-05 

44.6242981 

0 

0 

0 

89.2486057 
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TABLE II. - Continued. SAMPLE COMPUTATION OUTPUTS FOR PRINCIPAL MOMENTS OF INERTIA 

(b) Type 2 - PC1 3 


PYRAMIDAL 

WITH BOND 

LENGTH- 2.04300 , AND 

ANGLE- 100. 

10 


NO. 

ATOM 

KISLI 

MOLECULAR WEIGHT 

UK JCP V 18 PI 
X 

109 AUG 1950. 
Y 

l 

I 

P 

30.9749999 

0 

0 

0 

2 

CL 

35.4569998 

1.56617413 

0 

1.31184886 

3 

CL 

35.4569998 

1.56617413 

0 

1.31184886 

4 

CL 

35.4569998 

0 

1.96533303 

-0.55795603 

THE GENERATED SYMMETRIC MATRIX IS 

I 

I 32.47*808 0 

0 

I 

I 

THE EIGENVALUES* LAMBDA* OF THE HATRIX ARE 
PRINCIPAL MOMENTS OF INERTIA (G-SQ CM) 

I 

I 

I 

I 

0 

44.488149 

12.625000 

I 

I 

I 

I 

LAMBDA ( 1 )• 1 • OE+39 - 32.4760076 
LAMBDA(2)»1.0E+39 - 32.4768052 
LAMBDA(3)*1. 0E+39 - 57.7581563 

I 

I 

0 

12.625000 

45.746817 

I 

I 

PROOUCT LAMBDAS* 1.06+1 17 - 60920.0078 


THE GENERATED MATRICES FROM THE JACOBI METHOD ARE AS FOLLOWS* 


THE DIAGONAL MATRIX (D1 THE EIGENVECTOR MATRIX (S) 


I 

I 32.4768076 

t 

0 

0 

I 

I 

I 

I 1.00000000 

0 

0 


I 0 

I 

32.4768052 

0 

I 

I 0 

0.72449513 

0.68927990 


I 0 
I 

0 

57.7581563 

I 

I 

I 0 
I 

-0.68927990 

0.72449513 



THE MATRIX (S)*lO)*IS TRANSPOSE), WHICH 
SHOULD EQUAL THE ORIGINAL MATRIX IS* 


32.4768076 

0 

0 

0 

44.4881449 

12.6249990 

0 

12.6249990 

45.7468128 


13 



TABLE II. - Continued. SAMPLE COMPUTATION OUTPUTS FOR PRINCIPAL MOMENTS OF INERTIA 

(c) Type 3 
(c-1) CC1 4 


TETRAHEDRAL-REGULAR 

WITH 

BONO LENGTH= 

.7600 BARTELL, 

JCP V23 P 1 854 OCT 

1955. 

NO. 

ATOM 

MOLECULAR WEIGHT 

X 

Y 

Z 

1 

C 


12.0109999 

0 

0.58666666 

0 

2 

CL 


35.4569998 

1.43703395 

0 

0.82967)94 

3 

CL 


35.4569998 

0 

0 

-1.65934388 

4 

CL 


35-4569998 

-1.43703395 

0 

0.82967194 

5 

CL 


35.4569998 

0 

2.34666663 

0 

THE GENERATED SYMMETRIC 
I 

I 48.625861 

I 

1 

MATRIX IS 

0 

0 

I 

I 

I 

THfc EIGENVALUES, LAMBDA, OF THE MATRIX ARE 
PRINCIPAL MOMENTS OF INERTIA (G-SO CM) 

LAMBDA! 1 )• 1 . OE + 39 = 4B. 6258612 

LAMBDA12)* 1.0E+39 = 4B. 6258607 

I 

I 

0 


48.625861 

0 

I 

I 

LAMBDA! 3 ) • 1 • OE *39 = 4B. 6258612 

I 

I 

0 


0 

48.62586 1 

I 

I 

PRODUCT LAMBDAS*!. OE + 1 1 7 = 114974.601 


THE GENERATED MATRICES FROM THE JACOBI METHOD ARE AS FOLLOWS, 


THE DIAGONAL 
48.6258612 

MATRIX ( D ) 
0 

0 

I 

I 

THE 

I 

I 

EIGENVECTOR MATRIX !S) 
1.00000000 0 

0 

0 

48.6258607 

0 

I 

I 

0 

1.00000000 

0 

0 

0 

48.6258612 

I 

l 

I 

I 

0 

0 

1.00000000 


THE MATRIX ($)*!D)*(S TRANSPOSE), WHICH 
SHOULD EQUAL THE ORIGINAL MATRIX IS, 

1*8.6258612 0 0 

0 48.6258607 0 


0 


0 


48.6258612 



TABLE II 


Continued. SAMPLE COMPUTATION OUTPUTS FOR PRINCIPAL MOMENTS OF INERTIA 


(c) Concluded. Type 3 
(c-2) P 4 


TETRAHEDRAL-REGULAR WITH BOND LENGTH- 1.3533 MAXWELL JCP V3 P699 NOV 1935. 


NO. 

1 

2 

3 

4 

5 


ATOM 

MOLECULAR WEIGHT 

X 

Y 

Z 


0 

0 

0.45111433 

0 

P 

30.9749999 

1.10499991 

0 

0.63797200 

P 

30.9749999 

0 

0 

-1.27594399 

P 

30.9749999 

-1.10499991 

0 

0.63797200 

P 

30.9749999 

0 

1.80445731 

0 


THE GENERATED SYMMETRIC MATRIX IS 


THE EIGENVALUES, LAMBDA, DF THE MATRIX ARfc THE 
PRINCIPAL MOMENTS OF INERTIA (G-SQ CM) 


25.116959 0 


0 


0 


25.116959 0 


LAMBDA ( 1 ) • 1 . OE + 39 
LAMBDA ( 2 1 • 1 . OE + 39 
LAMBOAI 31*1. OE+39 


25. 1 169593 
25. 1 169593 
25. 1 169596 


0 


0 


25. 1 16960 


PRODUCT LAMBDAS* I.OE+1 17 = 15845.3263 


THE GENERATED MATRICES FROM THE JACOBI METHOD ARE AS FOLLOWS, 

THE 0 1 AGONAL MATRIX tD) THE EIGENVECTOR MATRIX (S) 


I 25.1169593 

0 

0 

I 

I 

I 

I 1.00000000 

0 

0 

1 0 

25. 1 169593 

0 

I 

I 0 

1.00000000 

0 

I 0 

0 

25.1169596 

I 

I 

I 0 
1 

0 

1.00000000 


THE MATRIX <S)*(D1*(S TRANSPOSE), WHICH 
SHOULD EQUAL THE ORIGINAL MATRIX IS, 

25.1169593 0 0 

0 25.1169593 0 


0 


25. 1 169596 



TABLE II, 


Continued. SAMPLE COMPUTATION OUTPUTS POP PRINCIPAL MOMENTS OP INERTIA 


(d) Type 4 
(d-1) (NaOH) 2 


NON-PLANAR WITH 3 CHAIN ATOMS, 1 PLANAR BRANCHES, AND 2 OUT-OF-PLANE ATOMS. 





ESTIMATES 

FROM NAOH AND L IOH, IL I OH) 2. 


NO. CHAIN ATOM 

BOND LENGTH 

ANGLE 

MOLECULAR WEIGHT 

X 

Y 

1 

NA 

2. 2500 



22.9909999 

0 

0 

2 

0 

3.4472 

40.00 


16.0000000 

-1.72360000 

1.44627209 

3 

0 




16.0000000 

1.72359994 

1.44627209 

BRANCH 

NOU 1*. ATTACHED TO CHAIN ATDM ( 

3 

, NUMBER 2 



NO. 

ATOM 

BOND LENGTH 

ANGLE 





4 

NA 

2.2500 

40.00 


22.9909999 

0 

2.692544 1 8 

O-P-A 

BRANCH NO. 

1, ATTACHED 

TD CHAIN 

ATOM 

0 , NUMBER 

2 


NO. 

ATOM 

BOND LENGTH 

O-P-A 

ANGLE 




5 

H 

0.9600 

120.00 


1.00799999 

-2.20359993 

1.44627209 

O-P-A 

BRANCH NO. 

2, ATTACHED 

TO CHAIN 

ATOM 

0 , NUMBER 

3 


NO. 

ATOM 

BOND LENGTH 

O-P-A 

ANGLE 




6 

H 

0.9600 

120.00 


1.00799999 

2.20359987 

1.44627209 


0 

0 

0 


0 


0. B3 1 38438 


-0.B313843B 


THE GENERATED SYMMETRIC MATRIX IS THE EIGENVALUES, LAMBDA, OF THE MATRIX ARE THE 

PRINCIPAL MOMENTS OF INERTIA (G-SQ CM) 


16. 199666 

—0.4701 138E— 07 

0.6131919 I 

I 

LAMBDA I 1 >*1.0E+39 

= 16.1778021 



I 

LAMBDA(2)*1.0E*39 

« 17.6397927 

0.470 1138E-07 

17.639793 

-0 I 

I 

LAMBDA(3)*1.0E»39 

- 33.3906244 

0.6131919 

-0 

I 

33.376763 I 

PRODUCT LAMBDAS*!. 

QE+117 = 9531.06799 


THE GENERATED MATRICES FROM THE JACOBI METHOD ARE AS FOLLOWS, 


THE DIAGONAL MATRIX (D) THE EIGENVECTOR MATRIX IS) 


16.1778021 

0.72469895E-18 

I 

0.5551 1 151E-16 I 

1 0.99936502 

-0.3211 1128E-07 

0.35630208E-01 

0. 72 46 9895 E- 16 

17.6397927 

-0.277S5576E-16 I 

I 0.321 35320E-07 

0.99999999 

-0.10629120E-09 

0.5551 115 1 E- 16 

-0. 2775S576E-1 6 

33.3986244 I 

I 

I -0.35630208E-01 
I 

0. 125 12 I I86-0S 

0.99936503 

THE MATRIX 4S)«(0)»(S TRANSPOSE), WHICH 
SHOULD EQUAL THE ORIGINAL MATRIX IS, 




16.1996632 

-0.4701 1378E-07 

I 

0.61319189 I 

j 




•0.4701 137BE-07 

17.6397922 

0.47184478E-15 I 

j 




0.613191 80 

-0. 16653345E- 15 

33.3767610 I 






TABLE II. - Continued. SAMPLE COMPUTATION OUTPUTS FOR PRINCIPAL MOMENTS OF INERTIA 
(d) Continued. Type 4 
(d-2) S 2 C1 2 


NON-PLANAR WITH 3 CHAIN ATOMS. 0 PLANAR BRANCHES. ANO 1 OUT-OF-PLANE ATOMS. 

BOWEN, CHEM. SOC., SPECIAL PUBL. NO. 11, 1958. 


NO. 

CHAIN ATOM 

BDND length 

ANGLE 

MOLECULAR WEIGHT 

X 

Y 

1 

CL 

1.9900 


35.4569998 

0 

0 

2 

S 

2.0500 

104.00 

32.0660000 

0.48142453 

1.93088849 

3 

S 



32.0660000 

2.53142446 

1.93088849 


I 


0 

0 

0 


O-P-A BRANCH NO. 1, ATTACHED TO CHAIN ATOM S , NUMBER 3 

NO. ATOM BOND LENGTH O-P-A ANGLE 

4 CL 1.9900 104.00 35.4569998 3.01284897 


1.93088849 -1.93088849 


THE GENERATED SYMMETRIC MATRIX IS 

I 

I 32.370342 -17.122919 17.122920 

I 

I 

1 -17.122919 54.089292 5.7624590 

I 

I 

I 17. 122920 5.7624598 54.089291 


THE EIGENVALUES, LAMBDA, DF THE MATRIX ATE THE 
PRINCIPAL MOMENTS OF INERTIA IG-SQ CM1 


LAMBDA! 1 1*1. OE+39 = 14. 8526815 
LAMBDA!2)*1.0E+39 = 65.8444834 
LAMBDA! 31*1.0E*39 = 59.0517456 


PRODUCT LAMBDAS*!. OE+1 17 = 50533.0405 


THE GENERATED MATRICES FROM THE JACOBI METHOD ARE AS FOLLOWS, 



THE DIAGONAL 

MATRIX ID) 


THE 

EIGENVECTOR 

MATRIX (S) 



14.8526815 

0 

I 

0.82881 94 2E-12 I 

I 

I 

0.81022295 

-0.50612171 

-0.2499064 1 E-08 


0 

65.8444034 

I 

0 I 

I 

0.41445062 

0.57291412 

0.707 10676 


0.82B8 1942E- 12 

0 

59.8517456 I 

I 

I 

-0.41445065 

-0.57291416 

0.70710672 


THE MATRIX <S>*!D1*(S TRANSPOSE), WHICH 
SHOULD EQUAL THE ORIGINAL MATRIX IS, 


32.3703346 

-17.1229141 

17.1229153 


-17. 1229138 
54.0892787 
5.76245928 


17. 1229153 
5.76245952 


54.0892787 



TABLE II. - Continued. SAMPLE COMPUTATION OUTPUTS FOR PRINCIPAL MOMENTS OF INERTIA 

(d) Concluded. Type 4 
(d-3 ) SF 4 


NON— PLANAR WITH 3 CHAIN ATOMS, 1 PLANAR BRANCHES, AND 2 OUT-OF-PLANE ATOMS. 

OODD TRANS FAR SOC V52 P1052 1956. FICTITIOUS THIRD ATOM. 


NO. CHAIN ATOM 

60N0 LENGTH 

ANGLE 

MOLECULAR WEIGHT 

X 

Y 

I 

F 

1.5800 



19.0000000 

0 

0 

2 

S 

0.7900 

60.00 


32.0660000 

-0.79000001 

1.36832011 

3 





0 

-0. 14901 161E-07 

1.36832011 

BRANCH 

NO. 1, ATTACHED TO CHAIN 

ATOM 

S 

, NUMBER 2 



NO. 

ATOM 

BONO LENGTH 

ANGLE 





4 

F 

1 . 5600 

60.00 


19.0000000 

0 

2.73664021 

0— P— A 

BRANCH NO. 

1, ATTACHED TO 

CHAIN 

ATOM 

S , NUMBER 

2 


NO. 

ATOM 

BOND LENGTH 

O-P-A 

ANGLE 




5 

F 

1.5600 

1 10.00 


19.0000000 

-1.33039176 

1 .366320 1 1 


I 


0 

0 

0 


0 


1.48471433 


O-P-A 

BRANCH NO. 

2, ATTACHED 

TO CHAIN 

ATOM 

S , NUMBER 2 

NO. 

6 

ATOM 

F 

BOND LENGTH 
1.5800 

O-P-A 
1 10.00 

ANGLE 

19.0000000 -1.33039176 


1.36832011 -1.48471433 


THE GENERATED SYMMETRIC MATRIX IS 


THE EIGENVALUES, LAMBDA, JF THE MATRIX ARE THE 
PRINCIPAL MOMENTS OF INERTIA (G-SQ CM) 


25.719417 -0 


0 


19.540759 -0 


LAMBDA! 1)* 
LAMBDA! 21* 
LAMBDA 131* 


1.0E+39 = 
1.0E+39 = 
1 . OE *39 * 


25.7194171 

19.5487505 

17.4537139 


0 


17.453714 


PRODUCT LAMBDAS*! ,0E*1 17 = 8775.42480 


THE GENERATED MATRICES FROM THE JACOBI METHOD ARE AS FOLLOWS, 


THE 0 1 AGONAL MATRIX ID) 
>5.7194171 -0 

0 

I 

I 

THE 

I 

I 

EIGENVECTOR MATRIX IS) 
1.00000000 0 

0 

) 19.5487585 

-0 

1 

I 

0 

1.00000000 

0 

) -0 

17.4537139 

I 

I 

I 

I 

0 

0 

1.00000000 


THE MATRIX (S1*!D1*!S TRANSPOSE), WHICH 
SHOULD EQUAL THE ORIGINAL MATRIX IS, 

25.7194171 0 0 

0 19.5487585 0 


0 


0 


17.4537139 



TABLE II. - Continued. SAMPLE COMPUTATION OUTPUTS FOR PRINCIPAL MOMENTS OF INERTIA 

(e) Type 5 - SOClg 


COORDINATES READ 

IN — U+VCOS (ALPHA) 

PALMER JACS 

V60 P2360 1938. 

S-0 1 .45A, S-CL 2 . 07 A , O-S-CL 106, CL-S-CL 

14. 

U, V, ALF 
ATOM 

INPUT X 

U 

V ALPHA 

U 

Y 

V 

Z 

ALPHA U V 

ALPHA 

S 

-0. 

0. -0. 

-0. 

-0. 

-o. -0. -o. 

-0. 

0 

-0. 

-0.7330 -0. 

-0. 

1.2506 

-o. -0. -0. 

-0. 

CL 

-0. 

2.0700 57.0000 

-0. 

-0. 

-0. -0. 2.0700 

33.0000 

CL 

-0. 

2.0700 57.0000 

-0. 

-0. 

-0. -0. -2.0700 

33.0000 

NO. 

ATOM 

MOLECULAR WEIGHT 

X 

Y 

z 


1 

5 

32.0660000 

0 

-0 

-0 


2 

0 

16.0000000 

-0.73303599 

1.25059398 

-0 


3 

CL 

35-4569998 

1.12740281 

-0 

1.73604809 


4 

CL 

35.4569998 

1.12740281 

-0 

-1.73604809 


THE GENERATED SYMMETRIC MATRIX IS 

I 

I 39.079304 4.3422745 

I 

I 

I 4.3422745 45.306801 

I 

-0 

-0 

I 

I 

I 

I 

I 

I 

THE EIGENVALUES, LAMBDA, OF THE MATRIX 
PRINCIPAL MOMENTS OF INERTIA (G-SQ CM) 

LAMBDA! 1 )»1.0E+39 * 36.0664145 
LAMBDA! 2 )• 1 • OE * 39 * 47.5997691 
LAMBDA! 3)*1.0E*39 * 13.4991376 

ARE THE 

I 

I 

-0 

-0 

13.499130 

I 

1 

PRODUCT LAMBDAS* 1 . OE* 1 1 7 = 23608.7295 


THE GENERATED MATRICES FROM THE JACOBI METHOD ARE AS FOLLOWS, 



THE DIAGONAL 

MATRIX (D) 



THE 

EIGENVECTOR MATRIX 

(S) 



36.8664145 

0 

0 

I 

I 

I 

I 

0.89096754 

0.45406698 

0 


0 

47.5997691 

0 

I 


-0.45406698 

0.89096754 

0 


0 

0 

13.4991376 

I 

I 

I 

I 

0 

0 

1 .00000000 


THE MATRIX (S)*<0)»(S TRANSPOSE), WHICH 
SHOULD EQUAL THE ORIGINAL MATRIX IS, 


39.0793824 4.34227443 0 

4.34227455 45.3867989 0 


0 


13.499 1376 



TABLE II. - Concluded. SAMPLE COMPUTATION OUTPUTS FOR PRINCIPAL MOMENTS OF INERTIA 

(f) Type 6 - S a Cl 2 


COORDINATES READ IN — 

X, Y, Z 

BOWEN, CHEM. SOC. 

SPECIAL PUBL. 

NO. 11, 1958. 

NO. 

ATOM MOLECULAR WEIGHT 

X 

Y 

Z 

1 

S 

32.0660000 

0 

0 

0 

2 

s 

32.0660000 

2.04999998 

0 

0 

3 

CL 

35.4569998 

-0.48142456 

1.93088858 

0 

4 

CL 

35.4569998 

2.53142080 

0 

1.93088849 

THE GENERATED SYMMETRIC 
I 

I 32.370343 

I 

I 

I 17.122910 

I 

MATRIX IS 

17.122910 

54.089228 

-17. 122889 
5.7624599 

I 

I 

I 

I 

I 

I 

THE EIGENVALUES, LAMBDA, 3F THE MATRIX ARE 
PRINCIPAL MOMENTS OF INERTIA (G-SQ CM) 

LAMBDA! 1 ) • 1 . OE + 39 = 14.8526877 
LAM8DA(2J* t.OE+39 = 65-8444166 
LAMBDA! 31*1. OE+39 = 59.8516831 

I 

I 

-17. 122889 

5.7624599 

54.089231 

I 

I 

PRODUCT LAMBDAS*!. 0E*1 17 = 58532.9443 


THE GENERATED MATRICES FROM THE JACOBI METHOD ARE AS FOLLOWS, 


THE DIAGONAL MATRIX (0) THE EIGENVECTOR MATRIX (SI 


14.8526877 

-0 . 888 1 7842 E- 1 5 

I 

—0. 4664 928 8E-09 I 

I 

I 0.81022272 

0.58612205 

-0.49451341E-06 

-0 . 88 8 1 7842E- 1 5 

65.8444166 

0.4547481 IE-12 I 

I -0.41445107 

0.57291483 

0.70710592 

-0 . 4664 9288E— 09 

0.4547481 16-12 

59.8516831 I 

I 

l 0. 41445067 
1 

-0.57291310 

0.70710757 


THE MATRIX <S)*(D)*(S TRANSPOSE), WHICH 
SHOULD EQUAL THE ORIGINAL MATRIX IS, 


I 

I 32.3703361 
I 

I 17.1229053 
I 

I -17.1228848 
I 

ENO OF FILE TAPE A 2 


17.1229050 

54.0892148 

5.76245928 


-17. 1228845 
5.76245952 
54.0892172 


20 



ooooooooooooonoononnooooonoooonononoonoooooonooooooonooo 


TABLE III 


CODE LISTING 


MOMENTS OF INERTIA 

PROGRAM COMPUTES PRINCIPAL MOMENTS OF INERTIA FOR THE FOLLOWING TYPE 


OF SUBSTANCES 

NTYPE= 1 PLANAR EFN 100 
NTYPE=2 PYRAMIDAL EFN 200 
NTYPE=3 TETRAHEDRAL EFN 300 
NTYPE=4 NON-PLANAR WITH OUT-OF-PLANE ATOMS, (O-P-A) EFN 400 
NTYPE=5 READ DATA (U, V, AND ALPHA), FOR EACH X,Y,Z EFN 500 
NTYPE=6 READ COORDINATES OF EACH ATOM, X, Y, Z EFN 600 


PERMANENT INPUT 

CARD 1 CC 1-3 NO. OF ENTRIES IN MOLECULAR WEIGHT TABLES. INTEGER 

CC 4-11 EN, AVOGADRO NO. *1.0E-23 FLOAT 

CC 12-19, 20-27, ETC. MOLECULAR WEIGHTS OF ATOMS, 

IN THE ORDER TO BE LISTED ON 
NEXT CARD. OVERFLOW ONTO 
ADDITIONAL CARDS IF NECESSARY, 

STARTING IN CC 1 WITH 8CC/ENTRY , 

AND MAXIMUM OF 10 WEIGHTS / CARD. 

CARD 2 NAMES OF ATOMS FOR WHICH MOLECULAR WEIGHTS HAVE BEEN 

READ IN. 2 CC / NAME, LEFT ADJUSTED, MAXIMUM 40 ATOM NAMES / CARD. 
OVERFLOW ONTO ADDITIONAL CARDS IF NECESSARY. 


INPUT FORMAT FOR EACH SUBSTANCE. 

CARD 3 CC 1-78 COMMENTS CARD. ALPH 

CC 79-80 THE 2 CHARACTERS Cl. 


CARD 4 CC 1-3 
CC 4-6 
CC 7-30 
CC 31-78 


NUMBER OF ATOMS. INTEGER 

TYPE OF MOLECULE, NTYPE=l-6 AS LISTED ABOVE. INTEGR 
NAME OF SUBSTANCE, EXAMPLE, BE10H)2 . ALPH 

NAME OF EACH ATOM, IN DESIRED ORDER, 

2CC / ATOM NAME, LEFT ADJUSTED, 

MAXIMUM 24 ATOMS. EXAMPLE, H 0 BEO H . 


THE REMAINING INPUT DEPENDS ON THE TYPE. 
BOND LENGTHS - ANGSTROMS. 

ANGLES - DEGREES. 


=1 , PLANAR. 


CARD 5 CC 1-3 NUMBER OF CHAIN ATOMS = N. INTEGER 

CC 4-6 NUMBER OF BRANCHES =M. INTEGER 

CARD 6 CC 1-10, 1 1-20 , ETC.,(N-1> BOND LENGTHS, FLOAT 

( N— 2 ) ANGLES, FLOAT 


OVERFLOW ONTO NEXT CARD (6A) IF NECESSARY (MAX 8 WDS./CD) 
ALL BOND LENGTHS ARE LISTED FIRST, THEN ALL ANGLES. 

IF THERE ARE ANY PLANAR BRANCHES, 

CARD 7 CC 1-3 NO. OF CHAIN ATOM TO WHICH BRANCH IS ATTACHED. INT 


MOM 001 
MOM 002 
MOM 003 
MOM 004 
MOM 005 
MOM 006 
MOM 007 
MOM 008 
MOM 009 
MOM 010 
MOM 0 1 1 
MOM 012 
MOM 013 
MOM 014 
MOM 015 
MOM 016 
MOM 017 
MOM 018 
MOM 019 
MOM 020 
MOM 021 
MOM 022 
MOM 023 
MOM 024 
MOM 025 
MOM 026 
MOM 027 
MOM 028 
MOM 029 
MOM 030 
MOM 031 
MOM 032 
MOM 033 
MOM 034 
MOM 035 
MOM 036 
MOM 037 
MOM 038 
MOM 039 
MOM 040 
MOM 041 
MOM 042 
MOM 043 
MOM 044 
MOM 045 
MOM 046 
MOM 047 
MOM 048 
MOM 049 
MOM 050 
MOM 051 
MOM 052 
MOM 053 
MOM 054 
MOM 055 
MOM 056 



i ii mu i 


TABLE III. - Continued. CODE LISTING 


CC 4-6 NO. OF ATOMS ON THIS BRANCH <N). 

CARO 8 CC 1-10, 11-20,. ETC., N BONO LENGTHS, N ANGLES. 

REPEAT CAROS 7 ANO 8 FOR EACH PLANAR BRANCH. 


CNTYPE=2 , 
C 

C CARO 5 


CNTYPE=3 , 
C 


PYRAMIDAL. 


CC 1-10 
CC 1 1-20 


BOND LENGTH. 
ANGLE. 


TETRAHEDRAL. 


INTEGER MOM 057 
FLOAT MOM 058 
MOM 059 
MOM 060 
MOM 061 
MOM 062 
MOM 063 
MOM 064 
FLOAT MOM 065 
FLOAT MOM 066 
MOM 067 
MOM 068 
MOM 069 
MOM 070 


c 

CARO 

5 

CC 

1-10 BOND 

LENGTH. 

FLOAT 

MOM 

071 

c 







MOM 

072 

c 







MOM 

073 

CNTYPE=4 

, 

NON-PLANAR WITH 

OUT-OF-PLANE ATOMS, 

(O-P-A). 

MOM 

074 

C 







MOM 

075 

C 

CARD 

5 

CC 

1-3 NO. OF 

CHAIN ATOMS. 

INTEGER 

MOM 

076 

C 



CC 

4-6 NO. OF 

PLANAR BRANCHES. 

INTEGER 

MOM 

077 

C 



CC 

7-9 NO. OF 

O-P-A BRANCHES. 

INTEGER 

MOM 

078 

C 







MOM 

079 

C 

INSTRUCTIONS FOR CARDS 

6, 7, AND 8, SAME AS 

FOR PLANAR, NTYPE= 1 • 

MOM 

080 

C 







MOM 

081 

C 

CARD 

9 

CC 

1-3 NO. OF 

CHAIN ATOM TO WHICH 

O-P-A BRANCH INTEGER 

MOM 

082 

C 





IS ATTACHED. 

MOM 

083 

C 



CC 

11-20 BOND LENGTH. 

FLOAT 

MOM 

084 

C 



CC 

21-30 BOND ANGLE. 

FLOAT 

MOM 

085 

c 



CC 

31-40 SIGN 

( *1.0 , OR -1.0), 

TO INDICATE ABOVE FLOAT 

MOM 

086 

c 





OR 

BELOW THE PLANE. 

MOM 

087 


REPEAT CARD 9 FOR EACH O-P-A BRANCH. 


CNTYPE=5 

C 


READ DATA CU, V, ALPHA), IN THE FORM U+VCOS ( ALPHA > , 

FOR X, Y, Z, OF EACH ATOM. 


CC 1- 8 U , TO CALCULATE 
CC 9-16 V , TO CALCULATE 
CC 17-24 ALPHA, TO CALCULATE 
CC25-32 U , TO CALCULATE 
CC33-40 V , TO CALCULATE 
CC41-48 ALPHA, TO CALCULATE 
CC49-56 U , TO CALCULATE 
CCS 7— 64 V , TO CALCULATE 
CC65-72 ALPHA, TO CALCULATE 


X OF GIVEN ATOM. 
X OF GIVEN ATOM. 
X OF GIVEN ATOM. 

Y OF GIVEN ATOM. 

Y OF GIVEN ATOM. 

Y OF GIVEN ATOM. 
Z OF GIVEN ATOM. 
Z OF GIVEN ATOM. 
Z OF GIVEN ATOM. 


REPEAT CARD 5 FOR EACH ATOM OF SUBSTANCE, N ATOMS=N CARDS. 


CNTYPE=6, 

C CARD 5 CC 1-10, 11-20, ETC., X OF EACH ATOM IN 

C CARD 6 CC 1-10, 11-20, ETC., Y OF EACH ATOM IN 

C CARD 7 CC 1-10, 11-20, ETC., Z OF EACH ATOM IN 

COMMON A , S , SD 

EQUIVALENCE I CHECK , NCHECK 1 , (NCOM,COMENT{ 14) ) 
EQUIVALENCE CA,D) 


COORDINATES READ IN, X, Y, Z. 


X OF EACH ATOM IN ORDER LISTED. 
Y OF EACH ATOM IN ORDER LISTED. 
Z OF EACH ATOM IN ORDER LISTED. 


MOM 088 
MOM 089 
MOM 090 
MOM 091 
MOM 092 
MOM 093 
FLOAT MOM 094 
FLOAT MOM 095 
FLOAT MOM 096 
FLOAT MOM 097 
FLOAT MOM 098 
FLOAT MOM 099 
FLOAT MOM 100 
FLOAT MOM 101 
FLOAT MOM 102 
MOM 103 
MOM 104 
MOM 105 
MOM 106 
FLOAT MOM 107 
FLOAT MOM 108 
FLOAT MOM 109 
MOM no 
MOM 111 
MOM 112 
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Continued. CODE LISTING 


DIMENSION A(3,3),D(3,3),S<3,3),SD(3.3) 

DIMENSION SUBS ( 4 ) , COME NT ( 1 4 ) , MOL C 26 ) , X ( 26 1 , Y« 26 ) , Z ( 26 ) , 
126 1 » WGHT ( 105), NAME! 1 05 1 , BOND ( 25 1 , ANG( 24 I , V ( 3) , ANGBR ( 2 1 


TO NUMBER OF INPUT TAPE 
TO NUMBER OF OUTPUT TAPE 


SET I T 1 

SET IT2 

I T 1*7 
IT2=6 

RAD=0. 0174532925 
READ PERMANENT DATA. 

READ INPUT TAPE IT 1 , 1 002, NTABL , EN, ( WGHT(J) , J= 1 .NTABL) 

NTABL-NTABL+1 

WGHT -(NTABL 1*0. 

READ INPUT TAPE I T 1 , 1 006 ,( NAME { J ), J* 1 , NTABL > 

READ INPUT CARDS 3 AND 4* INPUT FOR EACH SUBSTANCE. 

900 READ INPUT TAPE I T 1 , 1 003 , < COMENT (J ) , J* 1 , 14 ) 

901 READ INPUT TAPE I T 1 , 1 000 , NATMS, NTYPE, ( SUBS ( J 1 , J= 1 , 4 ) , 

1 ( MOLT J 1 , 1 ,NATMS ) 

WRITE OUTPUT TAPE I T2 , 20 1 5 , ( SUBS ( J 1 , J* 1 .4 > 

S W I TCH=0 • 

GET MOLECULAR WEIGHT OF EACH ATOM 

DO 801 J-l, NATMS 

IF ( MOL < Jl-NAMOLD ) 802*803, 802 

803 EMN l J ) =WTOLD 
GO TO 40 

802 DO 804 K«l, NTABL 

IF ( MOL ( J )— NAME ( K ) 1804,805,804 

805 KK*K 

GO TO 806 

804 CONTINUE 
GO TO 831 

806 EMN ( J l-WGHT (KK ) 

NAMOLD=NAME 1 KK 1 
WTOLD-WGHT(KK) 

40 EMt J)*=EMN< Jl/EN 
801 CONTINUE 

TRANSFER ACCORDING TO NTYPE, TO INPUT REMAINING DATA. 

GO TO 1100,200,300,400,500,600) , NTYPE 
ERROR WRITE OUT, AND GO ON TO NEXT CASE. 

831 CONTINUE 

WRITE OUTPUT TAPE IT2, 20 16,M0L( J * 

810 WRITE OUTPUT TAPE I T2 , 2020 , ( SUBS I I 1 , I* 1 ,4 ) 

B CHECK-230 1 60606060 


MOM 

113 

MOM 

l 14 

MOM 

1 15 

MOM 

116 

MOM 

1 17 

MOM 

118 

MOM 

119 

MOM 

120 

MOM 

121 

MOM 

122 

MOM 

123 

MOM 

124 

MOM 

125 

MOM 

126 

MOM 

127 

MOM 

128 

MOM 

129 

MOM 

130 

MOM 

131 

MOM 

132 

MOM 

133 

MOM 

134 

MOM 

135 

MOM 

136 

MOM 

137 

MOM 

138 

MOM 

139 

MOM 

140 

MOM 

141 

MOM 

142 

MOM 

143 

MOM 

144 

MOM 

145 

MOM 

146 

MOM 

147 

MOM 

148 

MOM 

149 

MOM 

150 

MOM 

151 

MOM 

152 

MOM 

153 

MOM 

154 

MOM 

155 

MOM 

156 

MOM 

157 

MOM 

158 

MOM 

159 

MOM 

160 

MOM 

161 

MOM 

162 

MOM 

163 

MOM 

164 

MOM 

165 

MOM 

166 

MOM 

167 

MOM 

168 
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CODE LISTING 


82 b RE AO INPUT TAPE I T 1 , 1 003 , ( COMENT ( J ) , J= I , 1 4 > 

DIFC=NCHECK-NCOM 

IF (DIFC*777777777777)825,901 ,825 
PLANAR SUBSTANCES INPUT 

100 READ INPUT TAPE I T 1 , l 005 , NCHA I N , NBR AN 

WRITE OUTPUT TAPE I T2 , 2 1 00 ,NCHA I N , NBRAN , ( COMENT C J ) , J= 1 , 13 ) 

105 DO 101 J= 1 , NA TM S 

101 Z ( J ) =0 . 

X { 1 ) =0 . 

v( n=o. 

N= 1 

IF ( NCHA I N— 2 ) 132, 102, 102 
132 WRITE OUTPUT TAPE I T2 , 20 1 8 , NCHA I N 
GO TO 810 

102 NMl=NChAIN-1 
NM2=NCHAIN-2 

READ INPUT TAPE I T 1 , 1 00 1 , < BOND ( J ) , J= 1 , NM U , { ANG ( J ) , J= 1 , NM2 > 

IF { ANG( 1) ) 115, 1 16, lib 

115 ANCLE=( 180.0-ANGt 1 ) )*RAD 

116 X ( 2 ) = { BOND ( 1) )*COSF (ANGLE) 

Y ( 2 ) =BOND ( 1 )*S INF (ANGLE) 

ANGDEG=0 • 

DO 1 12 J=3,NCHAIN 
ANGLE=ANGDEG*KAU 
SlGN=(-1 ) * • ( J- 1 ) 

X( J)=X ( J-l )+SlGN*BOND( J-l ) *C0SF( ANGLE ) 

IF (M0DF( ANGDEG, 1 80 . ) ) 1 1 4 , 1 1 3 , 11 4 

113 S I GN = 0 • 

114 Y( J)=Y( J-l )+SIGN*BOND( J-l ) *S I NF ( ANGLE ) 

112 ANGOEG=ANGDEG+ANG ( J- 1 ) 

WRITE UUTPUT TAPE I T2 , 2 1 08 , N • MOL ( 1 ) , EMN (1 ) , X ( 1 ) , Y ( 1 ) , Z ( 1 ) , BUND ( 1 ) 
IF (NCHA IN- 2 ) 1044,1044, 1045 

1045 WRITE UUTPUT TAPE I T2, 2 1 0 1 , ( J , MOL ( J ) , ANG ( J- l ) , L-MN ( J ) , X ( J ) , Y ( J ) , 
1Z( J),BGND( J),J=2,NMl) 

1044 N-NCHA I N 

WRITE CUTPUT TAPE I T2 , 2 1 02 , N , MOL ( N ) , EMN ( N ) , X ( N ) , Y ( N ) , l ( N > 

N = N+ 1 


MON 169 
MOM 170 
MOM 1 7 1 
MOM l 72 
MOM 173 
MOM 1 74 
MOM 17b 
MOM 176 
MOM 177 
MOM 178 
MUM 179 
MOM 1 BO 
MOM 181 
MOM 182 
MOM 183 
MOM 184 
MOM 1 85 
MOM 1 86 
MOM 187 
MOM 188 
MOM 189 
MOM 190 
MOM 191 
MOM 192 
MOM 193 
MOM 194 
MOM 195 
MOM 196 
MOM 197 
MOM 1 98 
MOM 199 
MOM 200 
MOM 201 
MOM 202 
MOM 203 
MOM 204 
MOM 205 
MOM 206 
MOM 207 


104 IF (NBRAN) 103, 135, 103 
103 CONTINUE 

DO 130 L - 1 , NBRAN 

READ INPUT TAPE I T 1 , l 005 , NATCH , NBRATS 
IF (NCHAIN- NATCH 1107,107,106 
107 WRITE OUTPUT TAPE I T2 , 202 1 , NC HA I N , N A TCH 
GO TO 810 

106 READ INPUT TAPE I T 1 , 1 00 1 , ( BOND ( J ) , J = 1 , NORA T S ) 
NP ! = N+ 1 
NN=NAT CH 
NNN= 1 
M=N 

128 ASSIGN 127 TO NNNN 

I F ( NATCH— 2 ) 121 , 122, 123 

121 ANGDEG= ANGBR ( 1 )-ANG( 1 ) 

GO TO 125 

122 ANGDEG = ANGBR t 1 ) 


MOM 208 
MOM 209 
MOM 210 
MOM 211 
MOM 212 
MOM 213 
MOM 214 

(ANGBR(J) ,J=1 ,NBRATS)MOM 215 
MOM 216 
MOM 217 
MOM 218 
MOM 219 
MOM 220 
MOM 221 
MOM 222 
MOM 223 
MOM 224 
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•b 


GO TO 125 

123 ANGOEG=ANGBR{ 1 > 

NAM 1=NATCH— 1 

DO 124 J-2 , NAM l 

124 ANGDEG=ANGDEG+ANG( J ) 

125 S I GN® I — 1*0) **N ATCH 

129 ANGLE=ANGDEG*RAD 

X { M )«=X ( NN ) +SI GN*BOND ( NNN ) »COSF ( ANGLE ) 

Y<M)=Y(NN)+SIGN*BOND<NNN)*S INF (ANGLE) 

IF I NBRATS-* 1 ) 145, 134,136 
136 GO TO NNNN, { 127, 134) 

127 ANGDEG=ANGDEG+ANGBR ( 2 ) 

SIGN-SIGN 

NN*N 

NNN=2 

M=M+1 

ASSIGN 134 TO NNNN 
GO TO 129 

134 WRITE OUTPUT TAPE IT2 , 2 1 03 , L , MOL ( NATCH ), NATCH 
DO 131 KK* 1 .NBRATS 

NN = N— 1 +KK 

131 WRITE OUTPUT TAPE I T2, 2 1 04 , NN, MOL { NN ) , BOND ( KK ) , ANGBRI KK ) , EMN ( NN ) , 
1 X ( NNs) , Y ( NN ) , Z ( NN ) 

130 N*N+NBRATS 

135 IFISWITCH1401»800,401 

145 WRITE OUTPUT TAPE I T2, 2 125 , NBRATS , L ,MOL< NATCH ), NATCH, MOL ( N ) , 
1MOHNP1 ) 

GO TO 810 

INPUT FOR PYRAMIDAL XY3 ALWAYS 4 ATOMS — X, Y, Y, Y 

BLN=X— Y BOND LENGTH, ALF* Y-X-Y ANGLE 

200 CONTINUE 

READ INPUT TAPE IT1 , 1001 , BLN, ALF 

WRITE OUTPUT TAPE IT2, 21 06, BLN, ALF, ( COME NT f J ) , J= 1 , 13) 

ANGL* ( ALF/2. ) *RAD 
S I NAL*S INF ( ANGL ) 

COSAL=COSF( ANGU 

SISQ*SINAL**2 

COSQ*COSAL**2 

xm-o. 

X { 2 )*BLN*S INAL 
X ( 3 )=— X ( 2 ) 

X ( 4 ) *0. 

Y< 1 ) *0 . 

Y(2>*0. 

Y ( 3 ) =0. 

TEM* ( 1*— 3«*SISQ+COSQ)/(2.*COSAL) 

YC4)=SQRTF(BLN**2*( 1.-TEM**2) ) 

U 11=0. 

1(2 )«=BLN*CCSAL 
Z(3)=Z<2) 

Z.( 4 ) = BLN*TEM 
GO TO 815 

INPUT FOR TETRAHEDRAL-REGULAR XY4 ALWAYS 5 ATOMS X Y Y Y Y 


MOM 225 
MOM 226 
MOM 227 
MOM 228 
MOM 229 
MOM 230 
MOM 231 
MOM 232 
MOM 233 
MOM 234 
MOM 235 
MOM 236 
MOM 237 
MOM 238 
MOM 239 
MOM 240 
MOM 241 
MOM 242 
MOM 243 
MOM 244 
MOM 245 
MOM 246 
MOM 247 
MOM 248 
MOM 249 
MOM 250 
MOM 251 
MOM 252 
MOM 253 
MOM 254 
MOM 255 
MOM 256 
MOM 257 
MOM 258 
MOM 259 
MOM 260 
MOM 261 
MOM 262 
MOM 263 
MOM 264 
MOM 265 
MOM 266 
MOM 267 
MOM 268 
MOM 269 
MOM 270 
MOM 271 
MOM 272 
MOM 273 
MOM 274 
MOM 275 
MOM 276 
MOM 277 
MOM 278 
MOM 279 
MOM 280 




1 
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C BLN=X-Y 80N0 LENGTH 
C 

300 CONTINUE 

READ INPUT TAPE IT1, 1001, BLN 

WRITE OUTPUT TAPE I T2* 2105 , BLN, { COMENT (J ) , J= 1 , 1 3 ) 
X( 1 )*0. 

X ( 2 )=BLN*0. 81 649658 
X ( 3 )*0» 

X ( 4 )=-X (2) 

X ( 5 ) =0* 

Y< I )=8LN/3* 

Y(2)«0. 

Y(3)*0. 

Y(4)=0. 

Y(5)=BLN*1. 33333333 
Z( 1)=0. 

Z(2)=BLN*0. 47140452 
Z { 3 )=— 2.*Z ( 2 ) 

Z ( 4 )=Z ( 2 ) 

Z ( 5 )*0 . 

GO TO 815 
C 

C NON— PL ANAR WITH OUT OF PLANE ATOMS 
C 


400 CONTINUE 

READ INPUT TAPE I T 1 , 1005*NCHAIN,NBRAN, NOPA 

WRITE OUTPUT TAPE I T2 , 2400 , NCHAI N, NBRAN, NOPA, ( COMENTI J ) , J= 1 , 13 ) 
$WITCH=0. 1 
GO TO 105 

401 DO 410 L* 1 » NOPA 

READ INPUT TAPE IT 1 , 1007* NATCH 1 8LN » ANGOEG* S IGNOP 

ANGLE-ANGDEG*RAD 

GO TO (402,402,403) *NATCH 

402 SIGN=+1.0 
GO TO 405 

403 S I GN=- 1*0 

405 X(N)«X (NATCH) +SIGN*BLN*COSF( ANGLE) 

Y ( N ) *Y ( NATCH) 

Z ( N )=S IGNOP*BLN*S INF ( ANGLE ) 

404 WRITE OUTPUT TAPE I T2, 2401 , L,HOL ( NATCH ), NATCH 

WRITE OUTPUT TAPE IT2*2 104,N,M0L I N) , BLN, ANGOEG,EMN( N) ,X (N) ,Y(N) , 

1 Z { N ) 

410 N=N+ 1 

GO TO 800 


C 

C COORDINATES READ IN — U+V*C0S ( ALPHA ) 

C 

500 CONTINUE 

WRITE OUTPUT TAPE IT2*2500» (COMENT ( J ) , J- 1 , 13 ) 

DO 501 L= 1 , NATMS 

READ INPUT TAPE I T 1 , 1500, ( BONDU ) , V ( J ) , ANG ( J ) , J= 1 , 3 ) 

WRITE OUTPUT TAPE I T2 , 250 1 , MOL (L ) , ( BOND ( J ) , V ( J ) , ANG < J ) , J= 1 , 3 ) 
X ( L ) =BGND ( 1 )+V( 1 )*(COSF(ANG( 1)*RAD) ) 

Y ( L )*B0ND { 2 )+V ( 2 ) * ( COSF ( ANG ( 2 ) *RAD ) ) 

Z ( L )=80ND( 3 ) +V ( 3 ) * ( COSF ( ANG < 3 ) *RAD ) ) 

501 CONTINUE 


MOM 281 
MOM 282 
MOM 283 
MOM 284 
MOM 285 
MOM 286 
MOM 287 
MOM 288 
MOM 289 
MOM 290 
MOM 291 
MOM 292 
MOM 293 
MOM 294 
MOM 295 
MOM 296 
MOM 297 
MOM 298 
MOM 299 
MOM 300 
MOM 301 
MOM 302 
MOM 303 
MOM 304 
MOM 305 
MOM 306 
MOM 307 
MOM 308 
MOM 309 
MOM 310 
MOM 311 
MOM 312 
MOM 313 
MOM 314 
MOM 315 
MOM 316 
MOM 317 
MOM 318 
MOM 319 
MOM 320 
MOM 321 
MOM 322 
MOM 323 
MOM 324 
MOM 325 
MOM 326 
MOM 327 
MOM 328 
MOM 329 
MOM 330 
MOM 331 
MOM 332 
MOM 333 
MOM 334 
MOM 335 
MOM 336 


I 


J 




f 
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GO TO 815 

MOM 

MOM 

337 

338 

COORDINATES READ IN X, Y, Z 

MOM 

339 

SPECIAL FOR MOLECULES NOT CATEGORIZED 

MOM 

MOM 

340 

341 

600 CONTINUE 

MOM 

342 

READ INPUT TAPE I T 1 , 1001 , ( X (J ) . J= 1 , NATMS ) 

MOM 

343 

READ INPUT TAPE IT 1 » 100 1 , ( Y ( J ) 9 J* 1 ( NATMS ) 

MOM 

344 

READ INPUT TAPE I T 1 , 1 00 1 , ( Z ( J ) , J= 1 , NATMS ) 

MOM 

345 

WRITE OUTPUT TAPE I T2 . 2600 , { COME NT { J ) , J= 1 , 1 3 ) 

MOM 

346 

815 WRITE OUTPUT TAPE I T2 *20 1 9 , ( J , MOL U ) , EMN I J } ,X I J ) t Y ( J ) , Z < J ) . 

MOM 

347 

1 J- 1 , NATMS ) 

MOM 

MOM 

348 

349 

COMPUTE ELEMENTS OF SYMMETRIC MATRIX, A. 

MOM 

MOM 

350 

351 

800 CONTINUE 

MOM 

352 

EMM=0. 

MOM 

353 

SMX=0. 

MOM 

354 

SMY=0. 

MOM 

355 

SMZ=*0* 

MOM 

356 

SXY=0. 

MOM 

357 

SXZ^O. 

MOM 

358 

SYZ=0. 

MOM 

359 

SMXY*0 • 

MOM 

360 

SM XZ = 0 • 

MOM 

361 

SM YZ=0 • 

MOM 

362 

DO 10 J= 1 , NATMS 

MOM 

363 

EMM=EMM+EM( J) 

MOM 

364 

EMX=EM( J)*X(J) 

MOM 

365 

EMY=EM( J)*YIJ) 

MOM 

366 

XSQ=X( J ) **2 

MOM 

367 

YSQ=Y ( J )#«2 

MOM 

368 

ZSQ=Z ( J ) 

MOM 

369 

SMX= SM X + EMX 

MOM 

370 

SMY=SM Y+EMY 

MOM 

371 

SMZ=SMZ+(EM(J)*Z(J) ) 

MOM 

372 

SMXY=SMXY+(EMX*Y( J) ) 

MOM 

373 

SMXZ~SMXZ+<EMX*Z< J) ) 

MOM 

374 

SMYZ=SMYZ+IEMY*Z( J) ) 

MOM 

375 

SX Y=SX Y + 1 EM ( J ) *( XSQ + YSQ) ) 

MOM 

376 

SXZ=SXZ+(EM(J)*(XSQ+ZSQ) ) 

MOM 

377 

SYZ=SYZ+IEM(J)*(YSQ+ZSQ) ) 

MOM 

378 

10 CONTINUE 

MOM 

379 

SMXSQ= SMX * *2 

MOM 

380 

SMYSQ= SMY**2 

MOM 

381 

SMZ SQ= SMZ* *2 

MOM 

382 

EY1 l=SYZ-( (SMYSO+SMZSQ ) /EMM ) 

MOM 

383 

EY22=SXZ-< (SMXSQ*SMZSQ)/EMM) 

MOM 

384 

EY33=SXY-( ( SMYSO+SMXSOJ/EMM) 

MOM 

385 

EY12 i =SMX*SMY/E MM—SMXY 

MOM 

386 

EY 1 3=SMX*SMZ/EMM— SMXZ 

MOM 

387 

EY23=SMY*SMZ/EMM-SMYZ 

MOM 

388 

A ( 1, U=EY1 1 

MOM 

389 

A{ 2 , 2 ) = EY22 

MOM 

390 

A ( 3, 3 ) =EY33 

MOM 

391 

A ( 1 , 2)*EY 1 2 

MOM 

392 


* 
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A ( 2 1 1 ) =EY 1 2 
A ( 1,31 =EY 1 3 
A ( 3 , 1 ) =EY 1 3 
A ( 2, 3 ) =EY23 
11 A ( 3 , 2 ) =EY23 


GET EIGENVALUES. 

CALL JACOBI (3, 1.0E-9) 

PROD 0=0 ( 1, l)*D(2,2)*DI3,3) 

WRITE OUT MATRICES AND EIGENVALUES. 


MOM 393 
MOM 394 
MOM 395 
MOM 396 
MOM 397 
MOM 398 
MOM 399 
MOM 400 
MOM 401 
MOM 402 
MOM 403 
MOM 404 
MOM 405 

WRITE OUTPUT TAPE I T2 , 2200 , E Y 11 , E Y 1 2 , E Y 1 3 , D { 1 , 1) , D ( 2 , 2 ) , E Y 1 2 , E Y22 , MOM 406 
1EY23,D(3,3),EY13, EY23 , EY33 , PRODD MOM 407 

820 WRITE OUTPUT TAPE I T2 , 220 1 , { ( ( D ( I , J ) , J= 1 , 3 > , ( S ( I , J ) , J= 1 , 3 )) , 1= 1 , 3 > MOM 408 
WRITE OUTPUT TAPE I T2, 2202 , < i SO ( I , J ) , J= 1 , 3 ) , I = 1 , 3 ) MOM 409 

MOM 410 

GO TO READ DATA FOR A NEW CASE. MOM 411 

MOM 412 

GO TO 900 MOM 413 

1000 FORMAT (213* 4A 6*25 A2) MOM 414 

1001 FORMAT ( 8F 1 0.4 ) MOM 415 

1002 FORMAT Il3»9F8.3/< 10F8.3)) MOM 416 

1003 FORMAT ( 1 3A6 * A2 ) MOM 417 

1005 FORMAT (313) MOM 418 

1006 FORMAT ( 40A2 ) MOM 419 

1007 FORMAT { I3,7X,3F10.4) MOM 420 

1500 FORMAT I 9F8 . 4 ) MOM 421 

2015 FORMAT ( 1H1 , 15X t 33HPRINCIPAL MOMENTS OF INERTIA OF 4A6 ) MOM 422 

2016 FORMAT ( 1 3HL THIS ATOM A6.36HIS NOT IN THE NAME AND WEIGHT TABLE. ) MOM 423 

2018 FORMAT { 52HL NO. OF CHAIN ATOMS IS LESS THAN 2. NCHA I N= I5JMOM 424 

2019 FORMAT(1HK,5X,3HNO.,6X,4HATOM,7X, 16HMOLECULAR WE I GHT , 6X , 1 HX , 1 6X , MOM 425 

1 1HY, 16X, 1HZ/I 1HJ, I7,8X,A6,4X,4G17.8) > MOM 426 

2020 FORMAT (30HL THIS CASE, SUBSTANCE IS 4A6,48H, CANNOT BE CONTINUMOM 427 

1 ED , GO ON TO NEXT SUBSTANCE. ) MOM 428 

2021 FORMAT ( 58HL BRANCH IS ATTACHED TO LAST OR NO CHAIN ATOM. THERE AMOM 429 
IRE I 3, 54H CHAIN ATOMS, THE BRANCH IS ATTACHED TO CHAIN ATOM N0.I31M0M 430 

2100 FORMAT ( 14HK PLANAR WITH I3,16H CHAIN ATOMS AND I3,10H BRANCHES., MOM 431 

14X, 13A6) MOM 432 

2101 FORMAT { I8*8X,A6, 15X,F9.2*4X,4G17.8/F35.4) MOM 433 

2102 FORMAT { I8,8X,A6,28X ,4G 17.8) MOM 434 

2103 FORMAT! 1HK,6X,1 OHB RANCH NO. I3.26H, ATTACHED TO CHAIN ATOM A6, MOM 435 

1 9H, NUMBER I 3/ 1 H J , 5X , 3HNO. ,6X , 4H ATOM , 7X , 11 HBOND LENGTH, 4X, MOM 436 

25HANGLE /) MOM 437 

2104 FORMAT (I8»8X,A6,F13.4,F1 1.2, 4X, 4G 17.8) MOM 438 

2105 FORMAT ( 40HL TETRAHEOR AL-REGULAR WITH BOND LENGTH- F8.4 ,3X, MOM 439 

113A61 MOM 440 

2106 FORMAT ( 30HL PYRAMIDAL WITH BOND LENGTH^ F10.5,13H , AND ANGLE= MOM 441 

IF 1 0.2/40X » 13A6 , A2 ) MOM 442 

2107 FORMAT ( 1HJ,40X, 13A6,A2) MOM 443 

2108 FORMAT ( 1 HK , 5X , 60HNO. CHAIN ATOM BOND LENGTH ANGLE MOLECMOM 444 

1ULAR WEIGHT, 7X, 1HX, 16X, 1HY, 16X, 1HZ/1HJ, 17, 8X, A6 , 28X , 4G 1 7. 8/F35 .4 ) MOM 445 

2125 FORMAT 1 6 1 HL NO. OF BRANCH ATOMS IS =, OR LESS THAN 0. NBRANCH ATMOM 446 
1 OM S= 1 4, 20H THIS IS BRANCH NO. I4/26H, ATTACHED TO CHAIN ATOM 6, MOM 447 
29H, NUMBER I4,30H. THE ATOMS ON THIS BRANCH ARE 2A6 ) MOM 448 
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2200 FORMAT ( 44HL THE GENERATED SYMMETRIC MATRIX IS , 36X , MOM 

2 46HTHE EIGENVALUES, LAMBDA, OF THE MATRIX ARE THE/80X , MOM 

3 39HPR INC I PAL MOMENTS OF INERTIA (G-SQ CM ) /8X , 1HI , 55X , 1 H I / MOM 

4 8X, 1HI,3G17*7,4X, 1HI/8X, 1 HI ,55X, INI , 18X, 1 9HLAMBDA (1) * 1 . OE + 39 = MOM 

5 G16.8/8X, IHI ,55X, IHI , 18X, 1 9HL AMBDAt 2 1 * 1 . OE+39 =G1 6.8/8X , IHI , MOM 

6 3G17.7,4X,1HI, 18X, 1 9HLAMBDA ( 3 ) * 1 . 0E + 39 =G 1 6* 8/8X , 1 HI , 55X , 1 H 1/ MOM 

7 8X , 1HI.55X, 1HI/8X, 1HI , 3G 1 7. 7, 4X , 1 H I , 1 8X, 26HPRODUCT L AMBD AS* 1 . 0E+ 1 MOM 

817 =G16.8/8X, 1HI,55X^1HI ) MOM 

2201 FORMAT ( 1H1 ,4X,61HTHE GENERATED MATRICES FROM THE JACOBI METHOD AREMOM 

1 AS FOLLOWS,/ 1HK,4X, 24HTHE DIAGONAL MATRIX (D) ,40X, 27HTHE E IGENMOM 
2VECT0R MATRIX (S) / 2H J 1 , 58X , 1 H 1 , 9X, 1 HI , 58X, IHI/ MOM 

3 2H 1,G16.8,2G20.8,2X,1HI, 9X, IHI ,G16.8, 2G20.8,2X, IHI/ MOM 

4 2H I ,58X, IHI , 9X, IHI ,58X, lHI/ MOM 

5 2H I,G16.8,2G20.8,2X, IHI, 9X, IHI , G 1 6. 8 , 2G20 • 8 , 2X , IHI/ MOM 

6 2H I , 58X , 1HI,9X,1HI , 58X, IHI/ MOM 

7 2H I ,G16.8,2G20.8,2X, IHI , 9X, IHI , G 1 6. 8, 2G20. 8 , 2X , IHI/ MOM 

8 2H 1 , 58X, IHI , 9X, 1 HI , 58X, IHI ) MOM 

2202 FORMAT ( 1HL , 4X , 39HTHE MATRIX (S)*(D)*IS TRANSPOSE), WHICH / MOM 

19X , 3 6H SHOULD EQUAL THE ORIGINAL MATRIX IS,/ MOM 

22HJI,58X, IHI/ MOM 

32H I,G16.8,2G20.8,2X, IHI/ MOM 

42H I ,58X , IHI/ MOM 

52H I,G 16.8,2G20.8,2X r IHI / MOM 

62H I., 58X, IHI/ MOM 

72H I,G 16.8,2G20.8,2X^1HI/ MOM 

82H I*S8X»1HI) MOM 

2400 FORMAT ( 1HL , 1 0X , 1 5HN0N- PLANAR WITH I3,14H CHAIN ATOMS, 13, MOM 

122H PLANAR BRANCHES, AND I3,21H OUT-OF-PLANE ATOMS. /1HJ,40X, MOM 

213A6J MOM 

2401 FORMAT ( 1 HK , 6X , 1 7H O-P-A BRANCH NO.I3.26H, ATTACHED TO CHAIN ATOM MOM 

1A6, 8H, NUMBER I 3/ 1 H J , 5X , 3HNO. , 6X , 4HAT OM , 7X , 1 1 HBOND LENGTH, 4X, 1 1HOMOM 
2-P-A ANGLE /) MOM 

2500 FORMAT l 1HL , 5X , 35HCOORD I N ATE S READ IN — U+VCOS I ALPHA ) f 5X , 1 3A6/ MOM 

1 1 HL , 1 9 HU , V , AL F INPUT X , 45X , 1 HY , 45X , 1 HZ /6H J ATOM, 6X , 1HU , 1 OX , 1HVMOM 

2^10X#5HALPHA, 1 6X , 1 HU * 1 OX , 1 HV , 1 OX , 5HALPHA , 16X, 1HU, 10X, lHV, I OX, MOM 

35HALPHA / ) MOM 

2501 FORMAT (2X,A6,3F1 1 . 4 , 1 OX , 3F 1 1 . 4 , 1 0X , 3F 1 1 .4 ) MOM 

2600 FORMAT ( 1HL , 5X , 35HCOORDI NATES READ IN — X, Y, l 13A6) MOM 

END ( 0 , 1,0, 1,0,1, 1,0,0, 1,0,0,0,0,01 


449 

450 

451 

452 

453 

454 

455 

456 

457 

458 

459 

460 

461 

462 

463 

464 

465 

466 

467 

468 

469 

470 

471 

472 

473 

474 

475 

476 

477 

478 

479 

480 

481 

482 

483 

484 

485 

486 
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SUBROUTINE JACOBI I NN*-FACTOR) 

SUBROUTINE JACOBI, COMPUTES ALL EIGENVALUES AND THE EIGENVECTOR 
MATRIX OF REAL SYMMETRIC MATRICES 
CALLING PROGRAM MUST, 

1 HAVE COMMON Y,S*SD 

2 HAVE DIMENSION Y (3*3 ) , S ( 3,3 > , SO { 3,3 ) 

ORIGINAL MATRIX AU,J> IN ARRAY Y. 

3 CALLING SEQUENCE , 

CALL JACOBI (N, FACTOR), WHERE N= ORDER OF MATRIX 

FACTOR^ MAGNITUDE OF FINAL THRESHOLD 
N IS AN INTEGER, FACTOR IS A FLOATING POINT NO* 

RETURN IS MADE TO CALLING PROGRAM WITH, 

1 DIAGONAL MATRIX {D>. IN ARRAY Y, THE DIAGONAL ELEMENTS ARE THE 

EIGENVALUES* 

THE ORIGINAL MATRIX, CA{I,J>), IS DESTROYED IN THE 
COMPUTATION PROCESS. 

2 EIGENVECTOR MATRIX IS) IN ARRAY S. 

3 MATRIX (S*D*S TRANSPOSE), WHICH SHOULD EQUAL THE ORIGINAL 
MATRIX, IN ARRAY SD. 

COMMON Y„S,SD 
EQUIVALENCE CY,D) 

DIMENSION Y(3,3),D{3f'3),S(3,3), SD (3,3) ,SDPR{3,3) 

1 N=NN 

2 S I GM A=N 

DO 200 K=1,N 
DO 200 L=1,N 
S0{K*L)=0. 

SDPRIK, L)=-0. 

200 S(K,L>=0. 

DO 201 K=1,N 

201 S(K,K)=1.0 

3 I.ND=0 
NM 1=N— 1 
K=2 
SUM=0. 

DO 300 I=1*NM1 
DO 301 J=K,N 
301 SUM^SUM+Yi I,JV**2 
K=K+l 

300 CONTINUE 

4 VI=SQRTF<2.*SUM) 

IF{VI)5,80,5 

5 V^VI 

VF=VI*FACTOR 

6 V— V/SIGMA 

7 J=2 

8 I*?l 

9 IF ( ABSFl YFI ,J ) )—V) 22^10, 10 

10 IND=l 

11 ALAM=-Y(I,J) 

12 AMU=I Y 1 1 , D— Y ( J , J ) ) /2. 


JAC 

001 

JAC 

002 

JAC 

003 

JAC 

004 

JAC 

005 

JAC 

006 

JAC 

007 

JAC 

008 

JAC 

009 

JAC 

010 

JAC 

on 

JAC 

012 

JAC 

013 

JAC 

014 

JAC 

015 

JAC 

016 

JAC 

017 

JAC 

018 

JAC 

019 

JAC 

020 

JAC 

021 

JAC 

022 

JAC 

023 

JAC 

024 

JAC 

025 

JAC 

026 

JAC 

027 

JAC 

028 

JAC 

029 

JAC 

030 

JAC 

031 

JAC 

032 

JAC 

033 

JAC 

034 

JAC 

035 

JAC 

036 

JAC 

037 

JAC 

038 

JAC 

039 

JAC 

040 

JAC 

041 

JAC 

042 

JAC 

043 

JAC 

044 

JAC 

045 

JAC 

046 

JAC 

G47 

JAC 

048 

JAC 

049 

JAC 

050 

JAC 

051 

JAC 

052 

JAC 

053 

JAC 

054 

JAC 

055 

JAC 

056 


TABLE III 


Continued 


CODE LISTING 


4 



IF (AMU ) 1 1 1, 112t 112 

JAC 

057 

1 1 1 

5GNMU=— 1 • 

JAC 

058 


GO TO 113 

JAC 

059 

112 

SGNMU=+1. 

JAC 

060 

113 

OMEGA=SGNMU*AL AM/ 1 SQRTF ( AL AM* *2+ AMU** 2 ) 1 

JAC 

061 


S I NTH=QMEGA/( SQRTF ( 2.* 1 l.+SGRTFl 1 .-0MEGA**2 ) ) ) ) 

JAC 

062 


SISQ=SINTH*SINTH 

JAC 

063 


COSTH= SQRTF (!• — SISQ) 

JAC 

064 


COSQ*COSTH*CO$TH 

JAC 

065 


SICO=SINTH*COSTH 

JAC 

066 


11=1 

JAC 

067 

13 

IFUI-I) 130,14 .130 

JAC 

068 

130 

I F C M-J)131,14. 131 

JAC 

069 

131 

B 1 P=Y( 1 1 ,.1 ) *COSTH- Y ( 1 1 , J)*SINTH 

JAC 

070 


B 1 Q=Y ( II, I.)*SINTH+Y(II ,J)*COSTH 

JAC 

071 


Y ( I l, I ) = 8 1 P 

JAC 

072 


Y 1 IL,J)=B1Q 

JAC 

073 

14 

SMIP=SII I , I ) *COSTH~S 1 I I, J)*SINTH 

JAC 

074 


SMIQ=S( I I, I )*SINTH+S( I I, J J*COSTH 

JAC 

075 


SC II . I ) = SM I P 

JAC 

076 


S( II..J ) = SM I Q 

JAC 

077 

15 

IFtN-lI ) 17, 17, 16 

JAC 

078 

16 

11=11+1 

JAC 

079 


GO TO 13 

JAC 

080 

17 

BP P=Y ( I , I ) +COSQ+Y ( J , J ) *S I SQ-2 • *Y ( I , J ) *S ICO 

JAC 

081 


BQQ=Y t I, I)*SISQ+Y( J, J 1 *COSQ+2 . *Y 1 I ,J)*SICO 

JAC 

082 


BPQ= ( Y ( I , I )-Y( J, J) 1 *S I CO+Y { I , J )* (COSQ-SISQ) 

JAC 

083 


Yl I, I ) =BPP 

JAC 

084 


Y ( J . J ) =BQQ 

JAC 

085 


Y 1 I , J ) =BPQ 

JAC 

086 

18 

I 1 = 1 

JAC 

087 

19 

Y ( I, LI ) = Y 1 II, I ) 

JAC 

088 


Y( J, II 1 = Y I II, J) 

JAC 

089 

20 

IFCN- 1.1)22, 22, 21 

JAC 

090 

21 

1 1 = I L+ 1 

JAC 

091 


GO TO 19 

JAC 

092 

22 

IF(I-( J-l) )23, 24,24 

JAC 

093 

23 

1=1 + 1 

JAC 

094 


GO TO 9 

JAC 

095 

24 

IFIN-J ) 26*26, 25 

JAC 

096 

25 

J=J+1. 

JAC 

097 


GO TO 8 

JAC 

098 

26 

IF { IND-1 )28,27,28 

JAC 

099 

27 

INO = 0 

JAC 

100 


GO TO 7 

JAC 

101 

28 

IF(V-VF)29,29,6 

JAC 

102 

IF OFF OLAGONAL ELEMENTS OF (D) MATRIX ARE LESS THAN 

1 • OE- 19, SET =0. JAC 

103 

29 

CONTINUE 

JAC 

104 


I-F4 ABSFCDl 1,2))- 1.0E-1 9) 950,951 ,951 

JAC 

105 

950 

Dt 1 , 2 ) =0« 

JAC 

106 


D ( 2, 1)=0. 

JAC 

107 

951 

IF ( ABSF ( D ( 1,3) J-l.OE- 19)952, 953, 953 

JAC 

108 

952 

D1 1 , 3 ) =0. 

JAC 

109 


013,11=0. 

JAC 

110 

953 

IF! ABSF 1012,3) )- 1 . OE- 1 9 ) 954 , 955 , 955 

JAC 

1 1 1 

954 

0(2,31=0. 

JAC 

1 12 


31 
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013,21=0. 

JAC 

1 13 

965 

CONTINUE 

JAC 

114 


DO 50 I* UN 

JAC 

115 


DO 50 J=UN 

JAC 

1 16 


DO 50 K=1,N 

JAC 

117 

50 

SDPRC I , J )=SDPR ( I , J ) +S ( I ,K ) *D( K, J ) 

JAC 

1 18 


DO 70 1= 1, N 

JAC 

119 


DO 70 J=UN 

JAC 

120 


DO 70 K= 1 , N 

JAC 

121 

70 

SD ( I , J USD ( I r J ) +SDPR ( I , K ) *S ( J,K ) 

JAC 

122 

71 

RETURN 

JAC 

123 

80 

DO 81 1=1, N 

JAC 

124 

81 

SDU,I )=Y t 1,1) 

JAC 

125 


GO TO 71 

ENDIO, 1, Or 1,0, 1,0,0,0,1,0,0.0,0,01 

JAC 

126 
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♦ DATA 

102 6.02322 


1.008 4.003 


6.94 


9.013 


10.82 


12.011 14.008 16.0 



19.0 

20.183 

22.991 

24.32 

26.98 

28.09 

30.975 

32.066 

35.457 

39.944 


39. 1 

40.08 

44.96 

47.9 

50.95 

52.01 

54.94 

55.85 

58.94 

58.71 


63.54 

65.38 

69.72 

72.6 

74.92 

78.96 

79.916 

83.8 

85.48 

87.63 


88.91 

91.22 

92.91 

95.95 

99.0 

101. 1 

102.91 

106.4 

107.88 

112.41 


114.82 

118.7 

121.76 

127.61 

126.91 

13L.3 

132.91 

137.36 

138.92 

140.13 


140.91 

144.27 

147.0 

150.35 

152.0 

157.26 

158.93 

162.51 

164.94 

167.27 


168.94 

173.04 

174.99 

178.5 

180.95 

183.86 

186.22 

190.2 

192.2 

195.09 


197.0 

200.61 

204.39 

207.21 

208.99 

210.0 

210.0 

222.0 

223.0 

226.0 

b. 

227.0 

232.0 

231.0 

238.0 

237.0 

242.0 

243.0 

247.0 

249.0 

251.0 


254.0 

253.0 

256.0 

254.0 








H HELIBEB C N 0 F NEN AMGALSIP S CLARK CASCTIV CRMNFECON ICUZNGAGEA5SE6RKRRBSR Y ZR 
N8MQTCRU RHP DAGCDJNSNSBTEI XECSBAL ACEPRNOPMSMEUGD T6DYH0ERTMYBLUHF TAW REOS I RPTAUHG 
TLPB8IP0ATKNFRKAACTHPAU NPPUAMCMBKCFESFMMDNQ 


WHITE JCP 
4 1H802 

4 0 


V 32 P488 FEB 1960. 


HOBO 


1.34 1.20 

JCP V32 P488 FEB. 
( HB02 ) 3 


1.36 


1.0 


1.0 


1.0 

JCP V 1 8 


1.36 


120 . 


120 . 


120 . 
1960 
8 0 

1.36 


240. 


1.0 
WHITE 
12 1 

6 3 

1.36 

120 . 

1 2 
i . 36 

3 2 
1.36 

5 2 

1.36 
K I SL I UK 

4 2PCL3 

2.043 100.1 

BARTELL, JCP V23 P1854 OCT 1955. 

5 3CCL4 C CLCLCLCL 

1.760 

MAXWELL JCP Vi P699 NUV 1935. 

5 3P4 P P P P 

1.353343 

ESTIMATES FROM NAUH AND L 1 OH , ( L I OH ) 2 . 

NAO 0 NAH H 


180. 


6-0 1.36A, 0-H l • A f CYCLIC 
6GBCJQHQHQH 


Cl 


Cl 


1.36 


120 . 


120 . 


120 . 


120 . 

Pi 109 


AUG 


240. 

240. 

1950. 

P CLCLCL 


Cl 


Cl 


Cl 


Ci 


6 
3 

2.25 
2 1 
2.25 
2 
3 

BOWEN 


4 { NAQH ) 2 
1 2 

3.4472 

1 

40.0 
0.96 
0.96 
CHLM. 


40.0 


4 

3 

1.99 

3 


4S2CL2 
0 1 

2.05 

1.99 


120.0 + 1.0 
120.0 - 1.0 
SOL., SPECIAL PUBL. 

CL S S 


NO. 

CL 


11, 1958. 


104. 

104.0 


- 1.0 


DOOO TRANS FAR SOC V52 P1052 1956. 


6 

3 

1.58 

2 

1.58 

2 

2 


4^F4 
1 2 


.79 

60. 
1.58 
i • :>8 


60. 


110 . 

110 . 


F S 


1. 

- 1 . 


FICTITIOUS THIRD ATOM. 
F F F 


Ci 


Cl 


i 


PALMER JACS V60 P2360 1938. 
4 5S0CL2 

0. 

-.733836 
2.07 57. 

2.07 57. 

BOWEN, CHLM. SOC., 

4 6S2CL2 
0.0 2.05 

0.0 0.0 1 , 

0.0 O.u 0- 


S-0 

S 


1.45A, 

0 CLCL 

L. 250594 


S-CL 2.07A, O-S-CL 106, CL-S-CL 114. 


2.07 

-2.07 


33. 

33. 


SPECIAL PUBL. NO. 

S S CLCL 

481424572. 5314208 
930888570.0 
0 1.9308885 


11, 1958. 


Cl 


Ci 


33 
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Figure 1. - Flow diagram of main program. 
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( C-1 ) Three-dimensional view. 
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Figure 3. - Concluded. Diagrams of certain molecules. 
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